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Abstract 

We include the effects of anisotropy and polarization in the hydrodynamics of inhomoge- 
neous vortex tangles, thus generalizing the well known Hall-Vinen-Bekarevich-Khalatnikov 
equations, which do not take them in consideration. These effects contribute to the mutual 
friction force F„s between normal and superfluid components and to the vortex tension 
force PsT. These equations are complemented by an evolution equation for the vortex line 
density, which takes into account these contributions. 

1 Introduction 

The possibility to use turbulent superfluids to explore with relative facility turbulent flows 
with high Reynolds numbers is an increasingly exciting perspective. Quantized vortices in 
superfluids have been mainly studied in two typical situations: rotating superfluids and coun- 
terflow experiments, the latter meaning the presence of a heat flux, but with zero bary centric 
motion. In these situations the vortices are modeled, respectively, as an array of parallel rec- 
tilinear vortices or as an almost isotropic tangle. In both cases, the mutual force between the 
normal component and the superfluid due to the presence of vortex lines is well known, and 
the so-called vortex line tension is zero 

However, in other situations, as rotating counterflow or non stationary Couette and Poiseuille 
flows, one expects a partially polarized vortex tangle, as a compromise between the orienting 
effect of a rotation or of a macroscopic velocity gradient, and the randomizing effect of the 
relative velocity of normal and superfluid components. In these cases, the mutual friction force 
Fns between these components, as well as the nonvanishing tension of the vortex lines PsT, 
which are fundamental ingredients of the hydrodynamics of turbulent superfluids in the well- 
known Hall-Vinen-Bekarevich-Khalatnikov (HVBK) model [5l[6], are not sufficiently known. 
Thus, the exploration of Fns and T for partially polarized tangles with nonvanishing average 
curvature of vortex lines is an open topic. 
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The quantitative values of the corrections are expected to be relevant in some steady 
states, where a sufficiently high polarization may induce macroscopic eddies [7], and many 
of the proposed experiments which involve nonstationary flows, for their own sake or in view 
to the application of classical measurement techniques, as for instance thin oscillating wires 
or small-particle velocimetry, to superfluids. Thus a basic understanding of these effects is 
essential to obtain the form of F^^ and T. This is the aim of this paper. 

The structure of the paper is the following one: in Section 2 an introduction to the usual 
HVBK model is made, with special emphasis on the expression of the mutual friction force 
and the tension T in the cases of rotating helium and counterflow superfluid turbulence; in 
Section 3 a microscopic expression for the friction force and vortex tension is given, checking 
their limit of validity in HVBK model; in Section 4 a more general expression of the HVBK 
model is given through more exhaustive expressions for the mutual friction force F^^ and 
tension T. The model is completed by a generalization of the Vinen's equation, as an evolution 
equation for the vortex line density L. At least, in Section 5 an application of the model to 
the interesting case of rotating counterflow turbulence is made, which takes advance of the 
recent study on the anisotropy of the vortex tangle. 

2 HVBK model. Evolution equations for the normal and the 
superfluid components 

It is known that if a superfluid (^He and ^He liquid helium, Bose-Einstein condensates, neutron 
stars, ...) rotates at a constant angular velocity exceeding a critical value Oc, an ordered 
array of quantized vortex lines of equal circulation n parallel to the rotation axis is created [1]- 
[4]. The quantum of vorticity is k = h/mi in superfluid ^He, with h the Planck constant and 
m4 the mass of ^He atom; in ^He and neutron stars, which are formed by fermionic particles, 
the superfluidity is due to the formation of Cooper pairs, and therefore it is k = h/2m, where 
m is the mass of an ^He atom or of a neutron, respectively. 

In most literature, the motion of a superfluid is modeled using Landau's two-fluid model, 
which regards the fluid to be made of two completely mixed components: the normal fluid and 
the superfluid, with densities p„ and ps respectively, and velocities v„ and respectively, with 
total mass density p and barycentric velocity v defined hy p = ps + Pn and pv = /O^Vs -|- Pn^n- 
The first component is related to thermally excited states (phonons and rotons) that form 
a classical Navier-Stokes viscous fluid. The second component is related to the quantum 
coherent ground state and it is an ideal fluid, which does not experience dissipation neither 
carries entropy. 

If the superfluid is put in rotation with angular velocity VL higher than Qc, the ordered 
array of parallel quantized vortex lines is described by introducing the line density L, defined 
as the average vortex line length per unit volume, equivalent to the areal density of vortex 
lines, which is proportional to the angular velocity [I]-[l], namely L = 2J7/k. 

It is well known too that a disordered tangle of quantized vortex lines is created in the 
so-called counterflow superfluid turbulence [I]-[lj, characterized by no matter flow but only 
heat transport, exceeding a critical heat flux Qc- When the turbulence is fully developed, the 
line density L is proportional to the square of the averaged counterflow velocity vector 
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related to heat flux (v„s = v„ — being the microscopic counterflow velocity) [T]-[l], L ~ 
llk^'ns / '^'^ ■> the dimensionless coefficient being dependent on the temperature. In (|2.ip 
and in the following, capital letters denote local macroscopic velocities averaged over a small 
mesoscopic volume A, threaded by a high density of vortex lines. 

Due to the smallness of the quantum of circulation, even a relatively weak rotation or a 
small counterflow velocity produce a large density of vortex filaments. It is therefore possible 
to develop a set of macroscopic hydrodynamic equations which average over the presence of 
many individual vortex lines and incorporate the macroscopic effects of the vortices in the 
evolution equations for superfluid and normal fluid velocities. 

A set of hydrodynamical equations frequently used is the Hall-Vinen-Bekarevich-Khalatnikov 
(HVBK) model [H El [6], [8]-[12]. These equations, which were derived by a number of re- 
searchers over the years, are those of the two-fluid model, modified to incorporate the presence 
of vortices. Here, for sake of simplicity, we consider the incompressible HVBK equations, 
which in an inertial frame are written [H [12] 

+ Pni^n ■ V)V„ = -^Vp„ - PsSVT + + T^V^V,, (2.2) 

ot p 

(9V n 

+ Ps(^^ ■ = --Vps + PsSVT - Fns + PsT. (2.3) 

ot p 

In these equations pn and ps are the effective pressures acting on the normal and the superfluid 

component, respectively, defined as Vp„ = Vp+ {ps/2)'S/V^^ and Vps = Vp — (/9„/2)Vl^^, p 

the total pressure, S the entropy, T the absolute temperature, and rj the dynamic viscosity 

of the normal component. The effects of the vortices are described by Fns, the friction force 

exerted by the superfluid component on the normal component — which bears an opposite 

sign in the evolution equation (j2.3p — and /OsT, the vortex tension force, related to the average 

curvature of vortices, which will be discussed below. 

The expression of the mutual friction force Fns in the HVBK equations is [12] 



FiHJBK) ^ p^^^ ^ [J ^ (y^^ -(3Vx u)] + p.a'Cd X (V,, - /3V x cj), (2.4) 

where a; = V x Vg is the local averaged superfluid vorticity, cj = the unit vector along 

uj, and P the vortex tension parameter, defined as [1] 



Kps 4:-K \aoL^/'^ 

with c a constant of the order of unity, ao the radius of the vortex core, of about 1 A, and ey 
the energy per unit length of vortex line. Coefficients a and a' depend on temperature and 
describe the interaction between the normal fluid and the vortices. They are linked to the 
well-known Hall-Vinen coefficients B and B' by the relations a = B{pn/2p), a' = B'(pn/2p) 
m. Note that this friction force is mediated through the presence of the vortex lines, as u) is 
related to them. 

In ()2.3p PsT is the vortex tension force whose microscopic meaning will be discussed in 
Section 3. In the HVBK equations T is substituted with 

rj^iHVBK) ^ (^y ^ ^) ^ ^ ^ . (2.6) 

From this expression it follows that in this approximation T vanishes when u) is homogeneous 
(Vlj = 0) or when the vortex lines are parallel to each other, in which case CJ will be orthogonal 
to the gradient of uj. 
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In the regular vortex array produced by the pure rotation, the vorticity uj equals 2ft 
everywhere, ft being the angular velocity, and the mutual friction force F„s (12. 4p assumes the 
Hall-Vinen expression 

= 2psan x[ftx (V„ - V,)] + 2psa'n x (V„ - V,), (2.7) 

because V x a; = 0, and uj = 2Ct. As seen in ()2.6p . the vortex tension T vanishes, because the 
vortices are straight lines. 

In counterflow superfluid turbulence the vortex tangle is supposed isotropic and the mutual 
friction force (j2.4p is expressed by [1]-[1] 

2 

Fns = --psHaLVns- (2.8) 

However, the hypothesis of complete isotropy of the tangle is not confirmed by the simulation 
of the dynamics of the vortex tangle made by Schwarz |13j-|15) and in the experiments |16j . 
Anyway, the assumed isotropy of the tangle implies that the vortex tension T vanishes because 
the microscopic vorticity would be isotropically distributed, in such a way that the average u) 
is null everywhere. 

The two quantities F^^ and T are the aim of the present paper. First, we will analyze the 
limit of validity of expressions (j2.4p and ()2.6p for Fns and T in the HVBK model, restricted 
to regions with a high density of vortex lines, all pointing in the same direction, and with the 
same curvature vector. Thus we will pay special attention to the fluctuations in the direction 
of the vortex lines and to the local anisotropy and polarization of the tangle as well as on its 
inhomogeneity, and their consequences on Fns and T. 

3 Expressions for F„s and T 

The HVBK equations ()2.2p ~ ()2.3p have been used in several problems in superfluid hydrody- 
namics, for instance, to study the Taylor-Couette flow in helium II |8]-|10t [T2 | and, recently, 
the jump in the rotational speed observed in neutron stars [17tll8j. The numerical simulations 
obtained with these equations show the formation of macroscopic vortices, but they do not 
include the microscopic chaotic dynamics expressed, for instance, by the fluctuations of the 
tangents to the vortex lines. 

Our aim here is to incorporate inhomogeneities, anisotropy, polarization, and tension of the 
tangle of vortex lines in an extension of the HVBK hydrodynamic equations. To this purpose, 
we make a critical analysis of these equations, and we perform a suitable modification of them. 
Furthermore, we include an additional equation for the evolution of L, besides the equations 
(j2.2p and (j2.3p . because L is related to forces and tensions. 

3.1 Microscopic expressions for friction force and vortex tension 

To derive the expressions (12. 4p and (12. 6p of Fns and of T, used in the HVBK equations, and to 
determine their limit of validity, we consider the microscopic description of a vortex tangle in 
the vortex filament model by Schwarz |13j-|15j. In this model, a quantized vortex filament is 
assumed as a classical vortex line in the superfluid with a hollow core and quantized circulation 
K. The vortex line is described by a vectorial function s(^, t), ^ being the arc-length, s' = ds/d^ 
the unit vector tangent along the vortex line, and s" = d'^s/d^'^ the curvature vector. The 
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normal component reacts to a moving vortex by producing the "microscopic" mutual friction 
force —(mf, which can be written as [l] 

(mf = apsKs' X [s' X (v„ - Vsi)] + a'psKs' X (v„ - Vsi), (3.1) 

where v^; = + Vi is the "local superfluid velocity", sum of the superfluid velocity at large 
distance from any vortex line and of the self-induced velocity, described by the tangent unit 
vector s' and by the curvature s". In the "local induction approximation", the self-induced 
velocity v; is approximated by [I]-[l] 

Vi ~ vi('°^) =p[s' X s"] . (3.2) 

The intensity of Vj is |vi| ~ P/R, with R the curvature radius of the vortex lines. The self- 
induced velocity is zero if the vortices are straight lines. Observing that Vn — v^; — ^ns — v;, 
being v„<j the microscopic counterflow velocity, the mutual friction force per unit length assumes 
the expression 

iMF = apsHs' X [s' X {Wns - Vi)] + a'psHs' X (v„s - Vi). (3.3) 

In the evolution equations for the velocities of normal and superfluid components, the 
"macroscopic" mutual friction force F„<j per unit volume which superfluid and normal compo- 
nents mutually exert, is obtained by averaging ()3.3p over a small volume A; one has 

< tMF >= r ,^ = TT / ^MFdt (3.4) 



fdC AL 

where the integral is made over all the vortices contained in the volume A; to obtain F„s, 
we must multiply the average ()3.4p . which denotes the average mutual friction force per unit 
length, by L, which denotes the length of vortex lines per unit of volume. One obtains 

Fns = [fns]av = L < ^Mp >= aps^L < s' X [s' X (v„s - Vi)] > +a' PsHL < s' X (Vns - Vi) > . 

(3.5) 

The vortex tension force T arises from the microscopic form of the evolution equation of 
Vg, which, neglecting the mutual friction force, is written [1] 

-Q^ = VL X LJmicr + V//, (3.6) 

where vl is the velocity of the vortex line and the term describes the effects of pressure 
and temperature gradients. Here denotes the microscopic velocity of the superfluid around 
the vortex and u^micr = V x v<j the vorticity of the single vortex line. 

Note that v^, = v^^ = + Vi [1] , with Vi the " self induced- velocity" . When the average 
on the mesoscopic volume A is taken one has 

— = [v, x (V X v,)],„ + [vi X (V X v,)]„, - + 5Vr, (3.7) 
ot p 

where Vp = —Vps/p + SVT. Furthermore, it is v.; x (V x v^) = — (v^ • V)vs + (1/2) Vv^. In 
the HVBK equation, one make the approximation [v^ x (V x Vs)]^^, = — (V^ • V)Vs, obtaining 

Ps-Q^ + Ps(Vs ■ V)V, = p,T - ^Vps + PsSWT. (3.8) 
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This allows to identify the tension T as 

T = [vi X (V X v,)]„„ = < Vi X s' >, (3.9) 

where we have denoted the average of vj x (V x v^), over the small vortex tangle contained 
in the mesoscopic volume A, with angular brackets, as in (13. 4p . 

Relation (13. 9p may be rewritten in several equivalent forms, by taking into account ()3.2p 
and some vectorial identities. For instance, < s' x Vi > may be expressed, in the local induction 
approximation, as 

< s' X Vi >= /5 < s' X (s' x s") >= -P < s" >, (3.10) 

where it has been taken into account that s' is a unit vector. 
Using ()3.10p . expression (j3.9p may be rewritten as 

T = [vi X OJmicAav = < (s' X s") X s' >= < s" > . (3.11) 

Another useful form for T may be obtained by using 

s" = (s' • V)s' = -s' X (V X s') (3.12) 

so that expression (13.101) may be rewritten as 

< s' X Vi >= /3 < s' X (s' X s") >= /? < s' X (V X s') > . (3.13) 

Expressions (13. 9p . ()3.1ip and (I3.13P will be used in the next sections. 

3.2 Limit of validity of the HVBK's equations 

As it can be seen by comparing (j2.4p and (j3.5p . in expression (j2.4p of the mutual friction force 
used in the HVBK equations the quantities < s' x (s' x (v„s — Vi)) > and < s' x [vns — Vi) > 
are approximated with Co x [u x {Yns — /3V x a))] and Co x (V„s — /9V x Co) respectively, and 
|V X Vs| is approximated by kL. The same approximations, when applied to the vortex tension 
T, yield (/3V xCb) xCj. 

We analyze in this section the limits of validity of the approximations leading to these 
expressions for the mutual friction force Fns and tension T, and we will show that these limits 
are related to the neglect of the second order moments of the vector s'. 

We consider first the two averaged quantities < s' x (s' x Vns) > and < s' x v„c, >. We 
denote with p =< s' > and V„s the averaged values of s' and v„s, and with Ss' and Svns their 
respective fluctuations. The average value of the unit vector s' is called tangle polarity [18j- 
[21j . and — recall that the superfluid vorticity is quantized — is linked to the local averaged 
superfluid vorticity Co by the relation 

For a totally polarized tangle with all the tangents s' parallel to each other, |p| = 1. In 
two dimensions, the polarization may be interpreted as (n"*" — n~)/(n~^ + n~) with n"*" the 
vortices rotating in one direction and n~ the vortices rotating in opposite direction. Processes 
producing a partial separation of + and - vortices will thus change p, and may lead to the 
formation of relatively large eddies [7]. Thus, the dynamics of p may be certainly important 
and it can be influenced by the boundary conditions. 
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With the notation just introduced, neglecting the fluctuations of the counterflow velocity 
V„s, we may write s' = p + Js'. Since < 6s' >= 0, we obtain 

< s' X (s' X v„s) >= p X (p X V„,)+ < 6s' X {6s' x V„^) >, (3.15) 

< s' X Vns >= p X V„s. (3.16) 

Note that in the equation (j2.4p only the first terms in the right-hand side of equations (j3.15p 
and (|3.16p appear. Indeed, recalling equation (j3.14p . the first terms in the right-hand side of 
equations (jS.lSp and (j3.16p can be written as 

p X (p X V„,) = i-ti X p X V„,] =< s' X (s' X v„,) yi^y^K)^ ^^ -^j^ 

P X V„, = \ujx Vns =< S' X V„, >(HVBK)^ (g -^g) 

kL 

whereas the other terms, quadratic in the fiuctuations, have been neglected. In general, this 
will not be correct, for instance, in the limiting situations of an isotropic tangle, p = 0, but 
< s' X (s' X v„,) >= (2/3) [19]. 

Note also that in this simplified hypothesis, the second term in equation (|3.15p . dependent 
on the second moments of s', can be neglected only if the orientational fluctuations of this unit 
vector in the small volume A are very small, i.e. if most of the vortex lines in the volume have 
the same direction. 

Furthermore, making use of equation (j3.12p . we obtain 

< Vi >= /3 < s' X s" >= /3 < (V X s') - [s' • (V X s')]s' >= /3 < (U - s's') • (V x s') >, (3.19) 

with U the unit matrix and s's' the diadic product. 

From equations (|3.17p ~ (j3.18p it is seen that in the HVBK equations all the second-order 
moments of the fluctuations in s' are neglected and it is implicitly assumed that kL ~ | V x | , 
that is I < s' > I ~ 1. Furthermore, we note that in the HVBK equations the quantity < Vi > 
is simply approximated by /3 < V x s' > and therefore the quantity < [s' • (V x s')]s' > is 
neglected. 

The relation 

kL ~ |V X Vsl = |a;|, (3.20) 

which is used to evaluate the line density L, is the most critical hypothesis underlying the 
HVBK equations. In fact, if in a mesoscopic region A there are several vortex lines oriented 
in a random way, the line density L in this point will be very different from L ~ |V x V^I/k, 
which corresponds to an extreme polarization with all or almost all the lines pointing out in 
the same direction, pointed out by | < s' > | = 1. For example, if in the small region A a 
vortex loop is present, the average circulation relative to this loop is zero, but it is not so for 
the line density L. 

As a consequence, the HVBK equations describe correctly the interaction between the 
normal component and the vortex tangle only in mesoscopic regions with a high density array of 
vortex lines, all pointing in the same direction (totally polarized) and with the same curvature 
vector. This fact explains the results of the numerical simulations obtained using the HVBK 
equations, which show the formation of macroscopic vortices and in which the microscopic 
chaotic dynamic behavior of the vortices does not appear. 



7 



4 Generalization of HVBK equations 



In a previous paper [SJj a hydrodynamical model of turbulent superfluids was formulated, 
which uses as fundamental fields the mass density p, the velocity v of the helium as a whole, 
the temperature T, the heat flux q and the line density L. In that work, only situations in 
which the tangle can be supposed approximately isotropic were considered, as is often made in 
the study of counterflow superfluid turbulence. In Ref. [19] the anisotropy of the vortex tangle 
was studied, restricting the study to homogeneous situations, and neglecting the influence of 
the vortex tension. In another paper p3] we have considered a more general situation, taking 
into account of inhomogeneities, anisotropy and polarization of the vortex tangle, studying 
the plane Couette and Poiseuille flow. For sake of simplicity, in that work, we have neglected 
the vortex tension T. 

4.1 New determination of Fns and T 

We want now to obtain expressions for F„<j and T overcoming some of the restrictions men- 
tioned in the Section 3.2. Using the local-induction approximation, and neglecting the fluctu- 
ations of the relative velocity V„,5, from (j3.5p we deduce that F„<j can be written |19j 



Fns = -PskL [a < U - s's' > +a' < W • s' >] • V^^ + ps^Lf] [a < s' x s" > +a' < s" >] , 

(4.1) 

with W the Ricci tensor (a completely antisymmetric third-order tensor such that W-s'-V„s = 
-s' X Yns)- 

Introducing the tensor H = + H", with [19] 



= ^ < U-sV >, n'^ = < W-s' >, (4.2) 
2 2 a 



the vectors I and J [H] 



(4.3) 



and ciL^/^ ~ XL I ^ characteristic measure of the vortex tangle introduced by Schwarz 

|15| . the mutual friction force (j4.ip can be written in a compact way as 



apsKL 



(4.4) 



The tensor H is especially useful to describe the geometrical properties of the tangle related 
with the orientational distribution of the vortex lines, whose local direction is indicated by the 
unit tangent s'. When the tangle is not completely polarized nor fully isotropic, it should be 
described by a tensor, rather than by a vector. The tensor U does not contain information on 
the curvature of the lines because it does not contain s"; this is furnished by the vectors I and 
J. 

The vortex tension T is also linked to the curvature vector J; indeed it results 

T = [vi X (V X vs)]av = kL< (/3s' X s") X s' >= kL/3 < s" >= kL^/^ciJ. (4.5) 

Observe that the tension of the vortex line, as the curvature vector, is zero in pure rotation 
with parallel straight vortex lines and in the isotropic tangle produced in well developed 
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counterflow turbulence, but it is not so in the presence of simultaneous counterflow and rotation 
[241 [25j or in the first stages of the turbulence |26l [27j or in the transient states after sudden 
acceleration in plane Couette and Poiseuille flow. For instance, T could be different than 
zero in Kelvin helical vortex waves, which appear when a sufficiently intense counterfiow is 
superposed to an axial rotation. In the transition from the array of straight lines to the array 
of helical vortex lines, a spiral tension would appear. The tension T would be a rotating 
quantity, as well as the curvature s" of the helical vortices. If the wavelength becomes shorter 
than the characteristic observational length, the average T will become equal to zero, but if 
the wavelength is high enough, it could lead to specific secondary flows of v^, through equation 

(1131). 

As a first modification to the HVBK equations we propose to take into account the fiuc- 
tuations of the vector s', which appear in the tensor < U — s's' >, but, for sake of simplicity, 
we propose to neglect the fiuctuations of s' x (V x s') = — s". 

To what concerns the vectors I and J, and the vortex tension T, with these approximations 
in mind, we have the following constitutive relations 



ci 



L^/^I =< s' X s" >~< U - s's' > -V X p, (4.6) 



ciL^/^j g// >^ -p X (V X p) = (p • V)p - ^Vp2. (4.7) 
The required equations for F„s and T are 

Fns = -PsKla < U - s's' > -[V^^ - /3(V x p)] + psKLa[p x Yns - (3px {V X p)] = 

= -^psfiLaU ■ [Vns - /3(V X p)], (4.8) 

PsT = -PskL(3p X (V X p). (4.9) 

In Section 5, some illustration of (j4.8p and (j4.9p will be presented, with different form 
of the field 11. But, to complete our hydrodynamical model, we must first add to equations 
j4.9p an evolution equation for the line density L. 



4.2 Generalized Vinen equation including polarization, anisotropy, and in- 
homogeneities 

In expressions ()4.ip and (|4.5p — or (|4.8p and (j4.9p — for F^,; and T, it appears L; therefore, to 
have a full description for the evolution of the system, an evolution equation for L is needed. 

The evolution equation for L in counterfiow superfiuid turbulence was formulated by Vinen 
[27j . Assuming homogeneous turbulence, such an equation is 

^ = ayVnsL^^^ - I3vkL^, (4.10) 

with a„ and (5^ dimensionless parameters. 

A microscopic derivation of this equation was given by Schwarz |13t [T5] on the basis of the 
dynamics of the vortices, neglecting however a term in the average curvature vector s". In 
Ref. [28J, because we were interested to study wall effects on the evolution of L, the following 
extension of equation ()4.10p was written 

^ ~ aciWns ■ 1^3/2 + q'ciV„, • - a/3c2L2, (4.11) 

at 
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with ci, I and J defined in equations (|4.3p . ci in the hne below (j4.3p . and C2 = / 

Substituting in this equation the relations (j4.6M.7l) . the following evolution equation for L 
which takes into account the polarization and the anisotropy of the tangle is obtained 

^ ~ aLVns- < U - s's' > -(V X p) - a'L\ns • p x (V x p) - ai3c2{p)L\ (4.12) 



In Refs. |23], [25]-[26], [29]-[3T] Vinen's equation was modified to describe more complex 
situations, as for instance the coupled situation of counterflow and rotation and in Couette 
and Poiseuille flows. Taking in mind the equation ()4.12p and the proposed equation of Ref. 
the previous equation would become 

^ 'l - ^ 



^ ~ aLVns- < U - sV > -(V xp)-a'LV„s-px (Vxp)-a/3KL2 
at 



1 - ^vIpI 

(4.13) 

where i? is a dimensionless coefficient lower than 1. A rigorous derivation of its form would 
require a microscopic extension of Schwarz's model including rotation effects. To check the 
consistency of the obtained evolution equation (I4.13|) for the vortex line density equation, in 
the next subsection an analysis based on the formalism of linear irreversible thermodynamics 
will be made. 

When inhomogeneities in the line density L are taken into account, the evolution equation 
for line density L must include a vortex density flux J-^, as [22] 



^ + V-3^ = a\ (4.14) 

where ai, stands for the right-hand side of equation (|4.13|) . The flux can also be expressed 
in terms of a convective part L\r^ with the tangle velocity and a dissipative part J^. The 
particular form of is also open to the debate. One can suppose that is an independent 
variable ^32J, or, more simply, supposing that is a dependent field: we have found [22J for it 
the form = uq^, vq being a coefficient describing the interaction between the vortex tangle 
and the heat flux q, which is linked to the counterflow velocity by the relation q = p^TSWus- 
The heat flux q may also be expressed in terms of VT and VL as [22j 

~- ''-VT-^VL. (4.15) 



Thus, the dissipative part of the vortex flux may be expressed as 

j^ = _^^Vr-^VL. (4.16) 
kLq kLo 

In isothermal situations, may be written as = —WVL, with D a vortex diffusion 
coefficient defined hy D = vqXq/i^Lq j28] . 



4.3 Consequences of Onsager-Casimir reciprocity relation 

In this Subsection, we show that another modification of the vortex line density evolution 
equation is necessary to insure the thermodynamic consistency of the evolution equations for 
L and for , according to the formalism of linear irreversible thermodynamics |31t [33t IM] . 
This analysis requires the presence of another term linked to the tension of vortices. 
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We follow the general lines of [HI [33l [34] with the aim to study the consequences of the 
Onsager-Casimir reciprocity relations on the evolution equations of and L proposed in 
the previous Sections. According to the formalism of nonequilibrium thermodynamics one 
may obtain evolution equations for and L by writing dVs/dt and dL/dt in terms of their 
conjugate thermodynamic forces —ps^ns and ey- The evolution equation for V^, neglecting 
inhomogeneous contributions of pressure, temperature and velocity, in an inertial frame, is 
written 

Ps—rr = -^ns - PsT = apsKL-n ■ V^s - apsnL$-Il • (V x p) + psKL(3p x (V x p) (4.17) 
dt 3 3 

and the evolution equation for L is written 

^ = aLVns- < U - s's' > -(V X p) - a'LV„, • p x (V x p) - q/3c2(p)L2. (4.18) 
dt 

However, in the right-hand side of M.lSh additional contributions must be included to make 
(j4.18p thermodynamically consistent with ()4.17p . 

Similarly to that presented in [31], we write dWg/dt and dL/dt in matrix form using 
the equations (|4.17p and M.lSp . and by means of Onsager-Casimir reciprocity we obtain an 
additional contribution to the evolution equation for L. The result is 

'ir\^r( -i^n -i:^^^ X p) . n + ^p X (V X p)\ /_^^v 

l^-^|(Vxp).n+^px(Vxp) -f^C2(p) ]\^v 

(4.19) 



where for C2(p) we can choose the expression C2(p) = ( 1 — \/|p| ) ( 1 — i^yTPl ) found in Ref. 
[23]. 

Therefore the equation for dL/dt becomes 



dL 

H 



aLYns- < U - s's' > -(V X p) - (1 + a')LV„, • p x (V x p) - a(3c2{p)L'^, (4.20) 



or, similarly, 

^ = ^aLVns • n • (V X p) - Xvns • T - «/3 (l - ^) (l - L\ (4.21) 

The new term not contained in the evolution equation ()4.18p for L is the coupling term 
between dL/dt and — /><jV„s in the matrix in (I4.19p . which is linked to the tension of vortices 
and it is null when the tension T is null. 

Note that, introducing the tensor 

1 -I- rv' 

n' =< U - s's' > < W • s' > (4.22) 

a 

the system (j4.19p can be written 

-Ifn -|;^n' . (V X p)\ /_ 

-|i^n'.(Vxp) -fLc2(p) I V 







(t) 


-A 




\ dt J 





(4.23) 
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5 Application to rotating counterflow turbulence 



Combination of counterflow and rotation is especially interesting in the context of the present 
paper, because it provides intermediate situations between an isotropic tangle and a totally 
anisotropic array of parallel vortices. The polarization and anisotropy of the tangle depend on 
the ordering influence of the rotation, which tend to align the vortices parallel to the rotation 
axis, and the randomizing aspects of the counterflow. In this section we will apply the previous 
results. 

5.1 Pure counterflow 

In most studies of counterflow superfluid turbulence, the vortex tangle is supposed homoge- 
neous and isotropic in the distribution of tangent vectors s' with respect to the counterflow 
velocity. Experimental observations and numerical simulations [15] show anisotropy in 
the vortex line distribution with vortices concentrated in planes orthogonal to V^^, and this 
anisotropy grows with the counterflow. Therefore, assuming V„s in the direction of the x axis 
and isotropy in planes orthogonal to it, one can choose for the tensor 11 (equation ()4.2|) ) the 
following expression [19] 

o / 2a \ 
= n|, = - 1 - a ; (5.1) 
^ V 1 -a / 

here a is the anisotropic parameter linked to the coefficients Jy and I± introduced by Schwarz 
|15j by the relations /y = 2o and I± = 1 — a. Being known the values of coefficient a, 
by accurate measurements of second sound, the second order moments of the unit vector s' 
remain determined. Indeed, it results 

< 42 >^ 1 _ 2a, < s'y^ >=< >= a. (5.2) 

For the vector I, we propose, according to (j4.6p . the following constitutive relation 

ciL^/^I =< s' X s" >~< U - sV > -V X p, (5.3) 

The asymmetric part of the tensor 11, the curvature vector J and the vortex tension T in 
this case are zero, owing to the supposed isotropy of the tangle, which implies < s" >= 0. 
Introducing (j5.ip into (|4.8p and (j4.13p we could obtain the influence of the anisotropy a on 
the mutual friction force Fns and the evolution of L. 

5.2 Simultaneous counterflow and rotation 

. Under the simultaneous influence of counterflow velocity V„s and rotation with angular 
speed ft, rotation tends to align vortex lines parallel to rotation axis, whereas counterflow 
velocity tends to produce a disordered (isotropic) tangle. In these situations one has partially 
polarized tangles, requiring the full detailed analysis presented here. We assume that the total 
ensemble of vortex lines is a superposition of both contributions [19] 

= (1 - b)njj + bn%, n" = cn]^. (5.4) 

In (15. 4p . b and c are parameters between and 1, describing the relative weight of the 
array of vortex lines parallel to and the disordered tangle of counterflow: when b = c = 
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we recover an approximately isotropic tangle, and when b = c = I the ordered array (total 
anisotropy). In general situations these coefficients depend on ft and V^s, or on the velocity 
gradient and V^s- 

In Ref. [E], supposing negligible the anisotropy due to the counterflow, as compared to 
the one produced by the angular velocity along the first axis, and assuming isotropy in the 
plane orthogonal to this axis, we obtained 

. - 6 

n^=| 1 + 1 |, n'^ = ^— c| 1 |, (5.5) 

1 + 1 




with the following interpretation of b and c in terms of moments of the tangent vector s' 

< 4 >= c, < 4 >= , < < >=< 4 >= ^— , (5.6) 



A microscopic evaluation of coefficients b and c was made in Ref. [19j, based on a para- 
magnetic analogy, which reflects the competition between the orienting effects of fi and the 
randomizing effects of V„s, respectively analogous to the orienting effects of a magnetic field 
H on magnetic dipoles /j, and the randomizing effects of thermal excitations. In rotating coun- 
terflow, the rotation 0, orients the vortices along its direction, in an analogous way to H, 
whereas the counterflow V^^ plays a disordering role. We found, using the Langevin model of 
paramagnetism, 

< 4 >= cothx , (5.7) 



and 

/2 2 ri 

< >= 1 H — cothx 

X 



(5.8) 



being V„ — along x axis. The first term on the right hand side, describes -j^^Fns, taking 



with X proportional to Qk/V^^. Similar situations would be found in Couette and Poiseuille 
flow where the velocity gradient, instead than a rotation, contributes to orient the vortices 
[551 ES] with Q replaced by the local shear rate. 

A particular illustration of the combination of ()5.5p and ()4.8p could be provided by the 
decay of small perturbations of V„ — V^. According to (j2.2p and (j2.3p and the assumption 
T = and VT = 0, one would have, in a linear approach 

7 = -\a^nL{l - 6)(V„ - V,) + ryV^V^, (5.9) 
at 3 Pn 

axis. 

into account (|4.8p and ()5.5p . 

For long-wave perturbations, V^V„ will be small as compared to the first term and the 
decay of V„ — to its steady state value will be exponential, with a relaxation time depending 
on the anisotropy parameter b. This time could provide a measurement of b independent of 
the measurement provided by the attenuation of the second sound along different axis of the 
system. An analogous analysis could be carried out to explore the anisotropy parameter a 
introduced in (15. ip for tangles in pure counterflow. 

A second particular illustration may underline the role of the term in I in expression (j4.4p 
for F„s. This would be reflected in the temperature gradient needed to maintain a steady 
state value of V„s, at constant pressure. Under these conditions, one has 

- 2psSVT + 2¥ns + ??V2V„ = 0. (5.10) 
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Thus, neglecting, for simplicity, the term in r/V^V„ and using (j4.4l) for F„<j, we have 



anL -11 • V„s + /3ciL 



^^/^I =SVT. 



(5.11) 



It is clear that the presence of I modifies the relation between V^^ and VT; this could allow 
to measure the influence of the corresponding term. 

These two simple illustrations show that the generalized equation for F„s considered in 
(14. 4|) is indeed expected to have specific applications. 

6 Conclusions 

Summarizing, in this work we propose to substitute in the HVBK equations the expression of 
the mutual friction force F„s (equation (|2.4p ) with the equation (j4.8p . to take into account of 
the second-order moment of s', allowing in this way that not all the vortex lines in the small 
volume element under consideration have the same direction. Analogously, the expression 
(12. 6p for the tension T has been replaced by equation (I4.9p . The coefficients b and c appearing 
in ()5.5p can be related to the counterflow velocity and to the angular velocity by the relations 
(15. 6p , with < > and < > expressed by (15. 7p and (j5.8p . For the vectors I and J we have 
chosen the constitutive relations (|4.6p and (j4.7[) . 

Another important modification consists in adding to the evolution equations for V„ and 
Vs an evolution equation for L including the effects of polarization anisotropy and inhomo- 
geneities. In fact, in a general situation, it is not correct in (I3.20|) to substitute kL with the 
modulus of the curl of V^. Relations (I4.13|) and ()4.14p avoid this simplification. 

In Section 5, we have concentrate our attention on rotating counterflow; another illustration 
of the possible physical consequences of the tension T, may be found in steady Poiseuille flows 
in an isothermal situation. There, addition of (j2.2p and (|2.3p leads for the equation describing 
the velocity profile of the normal component 



Thus, if T = 0, the V„ profile will have the typical parabolic form of Newtonian fluid. 
Modifications of fluids velocity profile have been studied by Godfrey and Barenghi [36j. From 
()6.ip . these modifications would be related to the form of T. A completely polarized array 
of parallel rectilinear vortices or a completely isotropic vortex tangle would have T = and 
would not modify the form of the V„ profile. However, net tension different from zero could 
arise at relatively low fluxes due to the influence of vortices pinned to the walls, which in the 
presence of the flow would have a curvature opposite to the velocity of the flow. In particular, 
they would have an influence on a small-amplitude oscillating flow along a cylinder with pinned 
vortices. 

Finally, to mention yet another situation where the effects of the polarization are impor- 
tant, and for which experimental data have been recently available, is the eddy formation in 
two-dimensional counterflow in the presence of a transverse cylinder [7] . The presence of the 
cylinder seems to separate clockwise and counterclockwise vortices, thus modifying the polar- 
ization. The non vanishing polarization turns out to be able to produce macroscopic eddies 
near the cylinder (downflow with respect to the velocity of the superfluid component), as a 
collective effect of the vortices. In this example, small changes in the polarization produce 
dramatic changes in the hydro dynamic flow. 



Vp + r/V^V^ + psT = 0. 
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